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ABSTRACT 

In this Letter I use the variation of the spread in rotation measure (RM) with Galactic latitude 
to separate the Galactic from the extragalactic contributions to RM. This is possible since the 
latter does not depend on Galactic latitude. As input data I use RMs from the catalogue by 
Taylor, Stil, and Sunstrum, supplemented with published values for the spread in RM ('ctrm') 
in specific regions on the sky. I test 4 models of the free electron column density (which I will 
abbreviate to 'DMo,') of the Milky Way, and the best model builds up DM„ on a characteristic 
scale of a few kpc from the Sun. <x RM correlates well with DM^. The measured cr RM can be 
modelled as a Galactic contribution, consisting of a term o~ rm,mw that is amplified at smaller 
Galactic latitudes as l/sin|fo|, in a similar way to DMoo, and an extragalactic contribution, 
0"rm,eg, that is independent of latitude. This model is sensitive to the relative magnitudes 
of ctrm MW and cr RM E g, and the best fit is produced by ctrm .mw ~ 8 rad/m 2 and <x RM ,eg ~ 
6 rad/m 2 . The 4 published values for ctrm as a function of latitude suggest an even larger 
0"rm,mw contribution and a smaller ct RM E g. This result from the NVSS RMs and published 
ctrm shows that the Galactic contribution dominates structure in RM on scales between about 
1° - 10° on the sky. I work out which factors contribute to the variation of <x RM with Galactic 
latitude, and show that the o~ R m,eg I derived is an upper limit. Furthermore, to explain the 
modelled cr RMMW requires that structure in (B\\) has a 1-cr spread < 0.4 ytzG 
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1 INTRODUCTION 

In recent years, models of the large-scale Galactic magnetic field 
have increased greatly in both complexity and breadth of input data 



and output variables ( s ee e.g. Fauvet et al l2010t Jaffe et all 



I i Nota & Katgertl |201CL liansson et aljl2009t IWaelkens et alj 12009 
Sun & Reichll2009t iMen et alj|2008l. [sun et alj|2008l. iBrown etal] 
2007). One popular approach to determine magnetic field strengths 
, is to look for Faraday rotation of t he polarized emissio n coming 
from extragalactic sources (see e.g. iGaensler et al]|2005h . The ro- 
tation measure (RM) quantifies the amount of Faraday rotation be- 
tween the source of the emission and us, the observers, and it de- 
pends on the free electron density, n c , and the magnetic field com- 
ponent along the line of sight, By, that the emission encounters 
along its path: 



RM [rad/m 2 ] 
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iree 



[cm' 3 ] By [jiG] dl [pc] 



where dl an infinitesimal part of the line of sight towards the ob- 
server. The RM can tell us about the magnetic field strength along 
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the line of sight, when the electron density contribution to RM is 
accounted for, for example by determining the free electron column 
density towards the source of the emission, which is known as the 
dispersion measure (DM) in the pulsar community. 



In this Letter I will focus on one particular observational as- 
pect, which is that the width of the RM distribution (cr RM ) of ex- 
tragalactic sources increases closer to the Galactic plane. This in- 
crease furthermore closely matches the increase in DM closer to 
the plane. I will use this, in combination with the fact that the ex- 
tragalactic RM contribution is independent of Galactic latitude, to 
separate this component from the Galactic RM contribution. This 
Letter is orga nised as follows. First I w i ll desc ribe the NVSS RM 
catalogue by iTavlor. Stil and Sunstruml ( 120091) that I use for my 
analysis in Sect. [2] and in Sect. [3] I describe the analysis that led 
to the observation that cr RM depends on Galactic latitude. In Sect. 
[4] I determine how well 4 models for the free electron density in 
the Galaxy predict the observed DM of pulsars at known distances. 
Finally, in Sect.[5]l will model the Galactic and extragalactic con- 
tributions to ctrm, and separate the two. Throughout this Letter I 
have calculated statistics in a way that is robust against outliers. 
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2 THE DATA 



iTavlor. Stil and Sunstruriil d2009h determined RMs for 37.543 
sources from the NVSS catalogue, which covers declinations > - 
40°. This means that the NVSS sources are separated on average 
by about 1° on the sky. Even though the accuracy of the individ- 
ual RMs is limited (the median error in RM is about 11 rad/m 2 ), 
the sheer size of the NVSS RM catalogue makes it very useful 
for studying the properties of the Galactic magnetized interstellar 
medium on large scales. In Fig.[T]l show the distribution of NVSS 
RMs as a function of Galactic latitude. 

Close to the Galactic plane, the RMs are much higher than fur- 
ther away, and this impacts the reliability of NVSS RMs at small 
Galactic latitudes. First, a large RM will induce a large amount of 
bandwidth depolarization, which reduces the polarized signal/noise 
ratio. Also, since the NVSS uses only 2 frequency bands, large RMs 
will suffer from nn ambiguities due to the periodicity of the polar- 
ization angles, and they will show up in the NVSS RMs as small 
RMs. Taylor et al. have suppressed the latter effect when they fitted 
RMs by also estimating the amount of bandwidth depolarization 
that would be expected for large RMs. Because of these reasons, 
I only use NVSS sources with |RM| < 300 rad/m 2 and Galactic 
latitudes \b\ > 20° in my analysis. The NVSS catalogue is not very 
sensitive to diffuse Galactic structure, which would have introduced 
further depolarization effects. 



3 THE SPREAD IN RM AS A FUNCTION OF LATITUDE 

Fig. shows that the spread in the NVSS RM increases for lines 
of sight closer to the Galactic plane, and in this Letter I will show 
that this increase very closely follows the increase in the Galactic 
DMoo . Fig. 5 from lTavlor. Stil and Sunstruml J2009h also shows this 
increase in o"rm- The NVSS RMs should however be corrected for 
2 effects before this conclusion can be drawn. 

First, the average RM will vary in a strip along Galactic longi- 
tude, an effect which becomes more pronounced close to the Galac- 
tic plane, and this increases <x RM when left uncorrected. I divided 
each strip along Galactic latitude into bins, and determined the av- 
erage RM for each bin. For latitudes within 77° of the plane bins 
are 5°/cos(fo) x 4° in size (Al x Ab). The l/cos(Z?) dependence of 
the bin width produces bins with a constant surface area on the sky, 
which guarantees that each bin will contain about 20 NVSS RMs, 
and that the maximum separation between two points in a bin is 
constant. Doubling Al does not significantly influence the results 
from Sect. [5] I divided the cap regions above \b\ = 11° into inter- 
vals with Al = 20° that converge at b=± 90°. I then calculated a 
cubic spline through the bin-averaged RM at the longitudes of the 
individual NVSS RMs, and subtracted this spline fit from the RMs. 
This removes structure in RM on scales larger than 10° (Nyquist 
sampling; for bins with \b\ < 11°). The resulting RM distribution 
for each strip along Galactic latitude is much better centred on 
rad/m 2 than the uncorrected measurements, and variations in (RM) 
that would increase cr RM are strongly reduced. 

Second, the ctrm that were determined from the corrected 
NVSS RMs should be corrected in a statistical sense for the mea- 
surement errors of the individual NVSS RMs. I estimate the cr RM 
that is produced by just the variation in uncertainty in RM of the 
NVSS sources ('errRM') by using a Monte Carlo simulation, which 
I set up as follows. First, I randomly draw errRM for 1000 lines of 
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Figure 1. Distribution of NVSS RMs as a function of Galactic latitude, 
which clearly shows the broadening of the RM distribution closer to the 
Galactic plane. The average RM and 1-cr spread around the average RM are 
calculated for 2° bins in Galactic latitude, and are shown as the red line and 
the blue lines on either side of the red line. Only lines of sight with |RM| < 
300 rad/m 2 are included in this and the following figures; 3% of the lines of 
sight have |RM| > 150 rad/m 2 and are not shown in this figure. 



sight in the NVSS catalogue, and I then use a random number gen- 
erato£] to find a distribution of RM based on these 1000 errRM. I 
then calculate the ctrm of this distribution, and repeat this process 
10 5 times. The total number of independently drawn random vari- 
ables for these 1000 lines of sight therefore does not exceed 10 8 , 
the maximum number of independent draws for the random num- 
ber generator that I used. I repeated this process for 35 sets of 1000 
lines of sight, to cover most of the errRM with \b\ > 5°. An analysis 
of the <r RM found after each run shows that the spread due to the 
errors in NVSS RMs gives a cr err RM = 10.4 ± 0.4 rad/m 2 . 

For 2 Gaussians, the widths (cr) add quadratically, and 



The randomn generator i n IDL, which is similar to the rani generator 
from Sect. 7.1 in lPress etal]jl992|) 
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Figure 2. Distribution of the excess ctrm as a function of Galactic lati- 
tude; circles (connected by a solid line) and squares (connected by a dashed 
line) indicate positive and negative Galactic latitudes respectively. The data 
points at Galactic latitudes of ±1° and ± 3° are missing, their ctrm values 
are 117 rad/m 2 and 85 rad/m 2 resp. The plotted data points contain at least 
15 RMs, and their RM distribution can be fitted by a Gaussian with a re- 
duced^ 2 < 2 (filled symbols) or between 2 <%L* < 4 (open symbols). The 
error bars (which are often smaller than the plot symbols) are calculated 
as the error in the mean RM of each latitude bin. Also shown are the un- 
corrected itrm calculated from Fig.[T](yellow line; sampled in 2° bins, and 
averaged over positive and negative latitudes) and the ctrm (not the excess 
o"rmO from Table [T] (red symbols). The red horizontal error bars indicate 
the range in Galactic latitude that is covered by these data sets. To illustrate 
how well the DM,*, distribution traces the excess ctrm at \b\ > 20°, I plot 
this distribution after scaling it by a factor of 0.37 (blue line; see Sect.EJ. 



°"errRM + nature = °"LrM + Mature- The ' excess <W. ^nature, Can then 

be recovered from the 

^"errRM + nature measured with NVSS, and 
the cr errRM that I simulated with a Monte Carlo process. A second 
Monte Carlo simulation showed that this relation is not exactly cor- 
rect due to a slight non-Gaussianity in the distribution of the RM 
errors in NVSS. This shows up as a 6% difference between the ex- 
pected o"rm excess and the measured ctrm excess, if the latter is 
10 rad/m 2 (this difference decreases with increasing cr rm excess). 
I corrected the cr rm excess for this effect, and I will use these cor- 
rected ctrm in the remainder of the Letter. 

In Fig. |2 1 plot the distribution of the excess <trm as a func- 
tion of Galactic latitude for both positive (solid line) and negative 
(dashed line) latitudes; in Fig. [3] I will distinguish between posi- 
tive and negative latitudes in the same way. I also show the ctrm 
distribution of the uncorrected NVSS RM data (yellow line). It 
is striking that there is a gap in ctrm between positive and nega- 
tive latitudes of about 2 rad/m 2 , but i ts origin is unkno wn to me. 
iTavlor. Stil and Sunstruml J200gh and IMao et all feOld) noticed a 
similar RM difference between the north and south Galactic caps. 

I compiled a list of published values of <trm in different re- 
gions on the sky, which I show in Table [T] and which I also in- 
clude in Figs. [2] — [3j These data sets are not limited by having 
only 2 frequency bands like the NVSS RMs. iFeain et alj J2009h 
removed a la rge-scale RM grad ient from their data before calculat- 
ing <x RM , and lMao et al.1 ( 120101) found an (almost) flat power spec- 
trum for their Galactic cap regions, so their ctrm is not apprecia- 
bly increased by a variation in RM. John ston-Hollitt et alj J2004I) 
argue that Galactic RM contributions are unimportant at latitudes 



region 


(b) 


o-rm [rad/m 2 ] 


source 


Centaurus A 


19° 


27±2 a 


Feain et al. (2009) 


LMC 


-33° 


15±l a 


Gaensler et al. (2005) b 


\b\ > 30° 


49 oc 


10±2 


Johnston-Hollitt et al. (2004) d 


\b\ > 77° 


81 oc 


9±0.3 a 


Mao et al. (2010) 



Table 1. Values of ctrm compiled from the literature. 

a : calculated as the standard error in the mean, using the number of sources 
that the authors listed to lie outside the target object. 
b : and Ann Mao, private communication 

c : (b) calculated as the surface-area-weighted average Galactic latitude of 
that data set 

d : and Melanie Johnston-Hollitt, private communication 

\b\ > 30°, and they would therefore not increase <trm. It turns out 
that there is some variation in Galactic RM even at these latitudes, 
which means that th eir ctrm is an upper limit to the excess ctrm. 
lGaensleretall ( l2005h did not correct for a smooth variation in RM, 
and also their ctrm is therefore an upper li mit. A variation in RM 
only i ncrea ses ctrm, but since the cr rm that Ijohnston-Hollitt et alj 
J2004) and lGaensleretalJ J2005I) calculate are very similar to the 
excess ctrm that I calculate from the NVSS RMs at these latitudes, 
their tr RM can be considered as tight upper limits. Since the data sets 
from Table Q] have (very) small measurement errors in RM, I will 
use t he ctrm, not the excess Crm, of these data sets. Since the data 
from I Johnston-Hollitt et al] d2004l) and the polar cap regions from 
IMao et al.1 feoifj) span a wide range in Galactic latitudes, I plot- 
ted these symbols at their surface-area-weighted average Galactic 
latitude^ which are (b) = 49° and (b) = 81° respectively. A uni- 
form source density will result in more sources at smaller Galactic 
latitudes, so using the surface area to weigh the Galactic latitudes 
reflects a weighting with the number of sources in each annulus. 

The correction for the variation in (RM) of the NVSS RMs is 
not perfect, and the RM distribution in the latitude bins often can- 
not be properly fitted by a Gaussian distribution, which is reflected 
in the high values for the reduced x 2 Ct^L) °f the Gaussian fits. 
However, Figs. [2] and [3] show that there is a reasonable agreement 
between the <Xr M from Table [T] and the excess cr rm that I calcu- 
lated, an indication that the corrections that I applied to the ctrm of 
the NVSS sources work well. Furthermore, latitude bins with x 2 ?d 
between 2-4 show the same global behaviour with latitude as the 
bins with xLa < 2. For my analysis in the remainder of this Letter 

1 will therefore use the excess cr rm of the NVSS sources with x 2 ed 
< 4, and the ctrm from the literature. I will distinguish excess cr RM 
from NVSS with x 2 cd < 2 from those with 2 < x 2 ei < 4 by using 
filled symbols for the former, and open symbols for the latter. 

4 A MODEL FOR THE FREE ELECTRON DENSITY 

Several models have been proposed in the literature to describe 
the free electron density in the Milky Way. Here I test the accu- 

2 I define the surface-area averaged Galactic latitude as 

90° 

f area in annulus X b db 

(b) = - — ^jjj for data with b > i m ; n . 

f area in annulus db 

<> > b mm 
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ra cy of the NE2001 mod el bv lCordes & Laziol feOO-tl 3 ! the model 
by_ Gaensler et alj J2008h . and the model by Berkhuijsen & Miillei 
(2008), by comparing the predicted DM to the measured DM for 65 
pulsars of which the distance is known to within 33%. I compiled 
this s ample from the ATNF pulsar catalogue (Man chester et al.l 
120051) . downloaded from this w ebsi ta__ on 25/05/2009, and the pa- 
pers bv lChatterieeetalJ d2009h and lDeller etaT] J2009I) . I also fit- 
ted a plane-parallel exponential model of the free electron density 
to this data set, taking into account the uncertainty in the pulsar 
distances. This model uses a more recent version of the catalogues 
that Gaensler et al. and Berkhuijsen et al. compiled, it includes also 
pulsars at latitudes within 40° from the Galactic plane (the model 
by Gaensler et al. is constrained by pulsars at latitudes \b\ > 40°), 
and it puts a stronger constraint on which pulsars to use (distance 
errors < 33%, where Berkhuijsen et al. required < 50%). The best 
fitting model has a mid-plane electron density /? c ,o = 0.02 ± 0.0001 
cm -3 and scale height h = 1.225 ± 0.007 kpc. (see also Schnitzeler 
& Katgert, in preparation) The ratios of the predicted DM/observed 
DM for these four models are 1.22 ± 0.53, 0.8 ±0.31, 1.05 ± 0.41 
and 0.99 ± 0.37 resp. (mean ± 1-cr spread; the error in the mean 
is smaller by a factor of V65 ~ 8), so I selected the final model 
to work with. The 65 pulsars that I selected cover a wide range of 
heights above the Galactic midplane, out to about 10 kpc, so the 
fitted scale height is well sampled by this pulsar distribution. 

The 3 models that use a plane-parallel exponential profile for 
the free electron density in the Milky Way (Gaensler et al., Berkhui- 
jsen et al., and the model I fitted) differ in mid-plane density and 
scale height. However, when integrating these models out to infinity 
(a good approximation for the NVSS sources), all 3 models should 
show the same asymptotic behaviour for a given Galactic latitude b: 
DM„ = 24.4/sin|£>|. DM(6=±90°) = 24.4 cnr 3 pc is the median DM 
that I calculated for the 8 pulsars fro m the sample that lie further 
than 4 kpc from the Galactic plane, teerkhuiisen & Mulled 120081 
find DM co (fe=±90°) = 21.7cirr 3 pc, and iGaensler et alj|2008l find 
DM„(£>=±90°) = 25.6 cirT 3 pc) In Fig._2]l overplot the DM„ pre- 
dicted by this model. To emphasize how well DM„ follows the 
trend in <x RM excess, I scaled the DM„ by a factor of 0.37, which 
is the average value of <trm X sin|£>|/24.4 for NVSS RMs at \b\ > 
20°. The resulting match in Fig. [2] is not perfect, since the values 
°f °"rm x sin|i| from Fig. [3] depend on latitude. In Sect. [5] I will 
show that the match can be improved by taking the fc-independent 
contribution from ctrm.eg to ctrm into account. 

To determine the characteristic distance over which most of 
the DM„ is built up, I introduce the electron-density weighted av- 
erage distance along the line of sight, (dist)„Jf| = h/s'm\b\. Note that 
(dist)„ c is somewhat larger than the distance over which half the 
DMo,, is built up, DM 1/2 =ln(2)x /i/sin|fc| ss 0.7 (dist>„ c . When going 
from |i|=90° to |i|=20°, (dist)„ e increases from 1 to 2h, or about 
1.2 - 2.5 kpc. These distances imply that the structure in DM„ 



3 and their website, http://rsd-www.nrl.navy.mil/7213/lazio/ne_model/ 

4 http://www.atnf.csiro.au/research/pulsaiVpsrcat/ 

5 I define the electron-density weighted average distance, (dist)„ e , as 

fine dl 

(dist)„ c = ^ , where both integrals are along the line of sight. For a free 

Jncd/ 



electron density that decreases exponentially away from the Galactic plane, 
<dist>„ c = h/sin\b\ 



(and tr rm) is neither produced very close to the sun, nor outside the 
Milky Way or its ionized halo. 

For a line of sight at 20° from the Galactic plane, (dist)„ e is 
2h, or about 2.5 kpc. This is nearby enough that the truncation of 
the Milky Way disk does not play a significant role in limiting the 
calculated DM„ . However, at smaller Galactic latitudes this can no 
longer be ruled out, at 10° from the plane (dist)„ c is already Ah (5 
kpc), and at 5°, (dist)„ c = 8/i (10 kpc), so the model I use to calculate 
DM,,, breaks down close to the Galactic plane. This might be the 
reason why in Fig.|2]DM„ increases much more rapidly at small 
latitudes than the measured <xrm. 



5 THE RELATIVE CONTRIBUTIONS OF THE MILKY 
WAY AND EXTRAGALACTIC SOURCES TO orm 

Fig. H] shows ctrm x sin|£>| as a function of Galactic latitude. The 
data points from Table [T] are nearly independent of latitude in this 
figure, which means that these cr RM closely follow a l/sin|£>| re- 
lation; the best-fitting DM„ model from Sect. [4] shows the same 
dependence on Galactic latitude. To explain why ctrm and DM^ 
might be correlated requires identifying which quantities determine 
how ctrm varies with Galactic latitude. 

The electron-density weighted line-of-sight component of the 
magnetic field, (B]|), is defined as 



(B\\) : 



j fi c S]| dl 

□nice 
observer 

J n e dl 



RM 



0.81 DM, 



The total RM of the line of sight can then be written in terms of its 
Galactic and extragalactic contributions as 

RM = RMmw + RM eg = 0.81<S||>DM„ + RM EG (1) 
with variance 

o-| M = 0.81 2 DM> 2 B||> + o.81 2 (fi„)V^ +o| MiEG (2) 

(assuming no correlation between the magnetic field and the free 
electron density) In this Letter <x RM is the (average) variance of 
an ensemble of lines of sight within a single cell used to remove 
the large-scale variation in (RM) from the NVSS data (Sect. [3) or 
within the area covered by the data sets from TableQ] 

A priori only <trm,eg is known to be independent of Galactic 
latitude; the first two terms on the right-hand side of Eqn. [2] can 
have both a ^-dependent and a ^-independent part. Although these 
different contributions to <x RM could in theory produce a compli- 
cated variation with latitude, the data show that this behaviour can 



be modelled in a simple way, as cr RM (b) ■■ 



li 'J'RM.MW \ 2 

y\ sin|i| ) 



.(I 



will use '<r' to distinguish the parameters from this model from the 
standard deviations in Eqn. [2} The model amplifies <trm,mw, which 
itself does not depend on Galactic latitude, at smaller Galactic lat- 
itudes as l/sin|£>| to mimic the increase in DM„ at these latitudes, 
and then combines it with (Trm.eg to form <trm. I show the best 
fitting models in Fig. [3] as blue lines, and these fits reproduce the 
global behaviour of <r RM with Galactic latitude (for \b\ > 20°). The 
parameters that produce these fits at positive (/negative) latitudes 
are cr RMMW = 6.8 ± 0.1 (8.4 ± 0.1) rad/m 2 and <Xrm EG = 6.5 ± 
0.1 (5.9 ± 0.2) rad/m 2 ; the values for the reduced^- 2 of these fits 
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Figure 3. Distribution of ctrm x sin|Z?| as a function of Galactic latitude, for 
the data points shown in Fig. [2] The data points at ± 1° and ± 3° are not 
shown; their ctrm x sin\b\ are 2.0 and 4.5 resp. I fitted the model described 
in Sect, \5\to the data points at \b\ > 20° for positive and negative latitudes 
separately, and show the best fits (solid blue lines). To illustrate the sensi- 
tivity of the modelled ctrm to <trm,mw , I also draw 2 curves on either side 
of the best-fitting curve to the data at negative latitudes, that only differ in 
0"RM,MW by 1 rad/m 2 from the best-fitting value, (dashed blue lines) 



are 4.4 and 3.9. I fitted the ctrm at positive and negative latitudes 
separately because they are off-set in Fig.[2] To illustrate how sensi- 
tive the modelled <r RM are to different (Trm mw, I show for negative 
latitudes also the best-fitting ctrm models for o"rm,mw=9.4 rad/m 2 
(top dashed curve) and o"rm MW =7.4 rad/m 2 (bottom dashed curve). 
Since the values of <t R m x sin|£>| depend much less on Galactic lat- 
itude for the data from table[T]than for the NVSS RMs, their ctrm 
can be modelled with a much smaller d~ rm.eg contribution. 

Now the (Trm.mw an d o"rm,eg parameters from the model 
can be identified with the terms in Eqn. [2] The DM„o"( B|| ) and 
(B^CTrjM^ terms can both have ^-independent parts, but since 
these terms add quadratically to <xrm,eg, o"rm,eg < ctrm,eg- The 
fc-dependent part of Eqn. [2] is modelled as <x RMMW /sin|£>|. When 
all the structure of the latitude-dependent part would be produced 
only by <x <fl|| ), which is equivalent to setting o"dm„=0, then a"< fi ) 
= o"R MMW /(0.81xDM oo (ii=90°)). cr DMoj will however not be zero, 
since it contains a contribution from the variation in modelled DM„ 
over the region of interest, and from local structure in the free 
electron density that is not included in the smooth DM„ model 
from Sect. [4] This structure in DM„ can explain part of the ob- 
served structure in ctrm, and the calculated <T( B|| > is therefore an 
upper limit, (again since the a terms add quadratically) Using 
DM 0O (£>=90°)=24.4 cirr 3 pc, and an average value of (Trm.mw = 7.6 
rad/m 2 , then gives o~( B|| ) ^ 0.4 \iQ. The <trm from TableQjgive the 
same value when calculating <t rm x sin|fc|/(0.81x24.4). (I ignored 
the 0"rm,eg term in this case - please see the end of the previous 
paragraph) 

The magnitude of the <trm,mw term compared to the ctrm,eg 
term implies that the Milky Way dominates structure in RM on 
scales between about 1-10 degrees, which can affect studies of 
the RMs of extragalactic sources, (the lower limit is set by the size 
scale probed by the NVSS RMs, and the upper limit is set by the 
size of the bins used to subtract the variation in (RM)) 



6 CONCLUSIONS 

In this Letter I show that the <t R m of NVSS sources can be mod- 
elled as a dominant Galactic contribution (cr RMMW « 8 rad/m 2 ) that 
is amplified at smaller Galactic latitudes as l/sin|£>|, similar to the 
increase in the Galactic free electron column density, DM„, and an 
extragalactic contribution that is independent of Galactic latitude 
(o"rm,eg ~ 6 rad/m 2 ). I corrected the ctrm of NVSS sources for the 
variation in (RM) along Galactic longitude, and for the broadening 
of the RM distribution that is produced purely by the uncertainties 
in the NVSS RMs. The 'excess' ctrm correlates well with DMc To 
calculate DM„, I compared 4 models of the free electron density, 
and I decided which one is the best based on how well it predicts 
the DM of pulsars at known distances. This model builds up DM„ 
on a characteristic scale of a few kpc for latitudes \b\ > 20°. I model 
the behaviour of <t rm with Galactic latitude as an extragalactic RM 
distribution, which does not depend on Galactic latitude, and the 
Galactic RM contribution, which is amplified at small Galactic lat- 
itudes to simulate the increase in DM„. The model follows the ob- 
served global behaviour in <t R m well, although there are significant 
localized deviations between the model and the observations. The 
best fitting values for the Galactic and extragalactic contributions to 
RM for positive (/negative) Galactic latitudes <xrm,mw = 6.8 ±0.1 
(8.4 ± 0.1) rad/m 2 and cr RM _ EG = 6.5 ± 0.1 (5.9 ± 0.2) rad/m 2 . By 
deriving which factors contribute to ctrm I show that this o"rm,eg 
is an upper limit. The spread in (B||> that is required to produce a 
°~rm,mw = 7.6 rad/m 2 is 0.4 fiG. However, fluctuations in the free 
electron density will also contribute to 0"rm,mw> which means that 
this <T( B|| > is an upper limit. I also compiled ctrm values from the lit- 
erature, that have more accurate RM determinations. These sources 
suggest that the o" RM MW contribution is even larger, and that o"rm,eg 
is very small by comparison. 

This result implies that structure in RM on angular scales be- 
tween about 1° - 10°is dominated by the Galactic foreground. 

Future radio polarimetry surveys will provide accurate rota- 
tion measures on a much finer grid on the sky. The analysis I present 
here can be improved upon with such data sets; in particular they 
will permit using smaller bins in Galactic longitude and latitude, 
which can better correct for the variation in RM with Galactic lon- 
gitude and latitude. 
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NOTE ADDED IN PRESS 

Since this paper was accepted it has become clear that the <x RM from 
Gaensler et al. and Johnston-Hollitt et al. can contain significant 
contributions from their errors in RM, which I should have taken 
into account. However, this does not affect the conclusions from 
this Letter. 
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